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THE APPLICATION OF HEAVISIDE’S STEP-FUNCTION 
TO BEAM PROBLEMS 


JOHN E. GOLDBERG, M. ASCE—Closure—The author wishes to thank the 
discussers for their critical and interesting comments and for their com- 
parisons of the step-function method with other and more general methods, 
including those founded upon the operational calculus. Admittedly, these 
techniques, such as the Laplace Transform, make possible the formal solu- 
tion of a much wider class of problems than the group treated by the author. 
Nevertheless, the essential simplicity of the step-function method, since it 
does not require the writing of a transform or an inverse, certainly has much 
in its favor for the type of beam problems treated in the paper. 

Complications arise in certain other types of problems, and for this rea- 
son the elementary step-function method was purposely applied only to a 
simple class of beam problems. The introduction of the sine transform of 
the step-function, as suggested by Professor van Langendonck, permits the 
relatively straightforward extension of the step-function technique to a wide 
class of beam-column problems and, consequently, should prove extremely 
useful in refined analysis. The use of the Laplace Transform, suggested by 
Mr. Saunders, instead of the elementary step-function alone, likewise is a 
much more general method, readily extended to buckling and beam column 
problems. It is, of course, true that, as stated by Mr. Saunders, rather 
comprehensive tables are available for finding the Laplace Transform of a 
great number of functions and the necessary inverse transforms. The author 
feels, however, that so long as the problem is of the type treated in the 
paper, the elementary step-function is sufficient to provide a solution in an 
extremely direct manner, and there is no need to employ the more general 
methods. Lest this statement be interpreted as unprogressive, the author 
hastens to state that he, too, probably would use a method related to the 
Laplace Transform for the more general problem. 

Professor Poggi has suggested and illustrated the use of the unit-impulse 
function in conjunction with the Laplace Transform, and has remarked upon 
the similarity of his method to that described by the author. In a sense, the 
unit-impulse function is a more flexible tool than the elementary step-func- 
tion and, consequently, might prove even more attractive to some engineers. 
Professor Poggi’s example, in which the Snitco-Hetenyi expansion is de- 
rived, is indeed interesting. One may remark, somewhat defensively, that 
expansion of the loading function, including concentrated loads, into a Fourier 
sine series is also extremely useful in many cases but, possibly, because of 
its transcendental nature, not quite as convenient as the Maclaurin series 
used by Professor Poggi. 

The author concurs with Professor van Langendonck regarding the re- 
striction upon the primitive, F(x), appearing in Equation (3). It does indeed 
seem necessary that F(O) = O for complete generality. This requirement is 
identically satisfied in the problems which appear in the paper. 
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DISCUSSION OF INELASTIC BEHAVIOUR OF REINFORCED CONCRETE 
MEMBERS SUBJECTED TO SHORT TIME STATIC LOADS 
PROCEEDINGS—SEPARATE NO, 286 


B. BRESLER, ! J.M., ASCE—Mr. Lee’s application of the inelastic flexure 
theory to reinforced concrete members is a valuable addition to the literature 
on this subject. An interesting feature of Mr. Lee’s paper is the assumption 
of parabolic stress-strain relationship for concrete in flexure and the verifi- 
cation of this assumption by beam tests. In deriving the relation between 
concrete stress and observed strains, Eq. (3), the author assumes an idealized 
stress-strain relationship for steel shown in Fig. 4 of the paper. For steel 
stresses between proportional limit and yield point this assumption may in- 
troduce significant errors if the actual steel stress strain diagram varies 
appreciably from that assumed, 

The importance of the actual shape of steel stress-strain diagram becomes 
apparent if values of concrete stress corresponding to steel yield strain es 
are computed from both Eqs. (3) and (4). Since data in Table 1 of author’s 
paper are limited to steel strains just below yield, the following estimated 
values will be used to illustrate the discrepancy between Eqs. (3) and (4). 
Consider yield strain e, equal to 1615 x 10%, and an estimated value of con- 


e 
crete strain e, equal to 2310x 10 s from which the value of Gen) is computed to be 


0.40. Using these values f, computed from Eq. (3) is 3600 psi, and that com- 
puted from Eq. (4) is 2130 psi. The discrepancy between these values of f, 
indicates that at strains approaching yield author’s Eqs. (3) and (4) are not 
valid. 

Author’s Eq. (2) can be rewritten in terms of steel stress f, as follows: 


(2') 


Substituting k in terms of e, and e, and differentiating both sides with respect 
to @¢, is: 


Ges ) ats 


+ (3) 


f. = pif, (it 


Equation (3') is the general expression for f{,, and author’s Eqs. (3) and (4) 
can be derived from it for the assumed shape of steel stress-strain diagram. 

The flexural stress-strain curve for concrete can be obtained from beam 
test data using either Eq. (2') or Eq. (3'). Differentiation of Eq. (2') with 
respect to concrete strain determines the value of extreme fiber stress f, 
corresponding to the extreme fiber strain e, as follows: 


1, Associate Prof. of Civ. Eng. Univ. of California, Berkeley, Calif. 
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The stress-strain relation for a given range of strains can be determined 
from Eq. (2") using either graphical or numerical differentiation. 

An experimental investigation of flexural stress-strain relationship based 
on Eq. (2") was carried out as part of M.S. Thesis? at the University of Cali- 
fornia, Berkeley. Three reinforced concrete beams 11 ft. 6 in, long having 
different concrete strengths, cross-sectional dimensions, and amounts of re- 
inforcement were subjected to third-point loading and tested to failure. 
Strains were measured at three levels in the concrete and also in the steel 
reinforcement, 

Stress-strain relationships for the three beams tested were determined 
from Eq. (2°) using numerical differentiation procedure. Three 6 x 12 in. 
cylinders were cast with each beam and these were tested to determine the 
stress-strain relation in compression, In Fig. 15 flexural stress-strain re- 
lationship obtained from one of the beams is compared with the average of 
three compressive stress-strain curves obtained from cylinder tests, It can 
be seen that the flexural stress-strain diagram varies considerably from the 
cylinder stress-strain curve in the region of high stress. 

It should be pointed out that calculation of flexural stress-strain relation- 


ship involves numerical differentiation which may be sensitive to the accuracy 


of experimental measurements, The experimental accuracy obtained in the 
tests described here is sufficient to allow the type of analysis proposed and 
therefore the curve shown in Fig. 15 is a reliable representation of test data 
obtained, Although test results described in this discussion seem to indicate 
that flexural stress-strain relationship for concrete may differ considerably 


from that obtained in compression, results of this limited investigation cannot 
be considered conclusive. Furthermore, the writer considers that the precise 
shape of the flexural stress-strain is primarily of academic interest as effects 


of plastic flow, shrinkage, size, previous history of stress and strain, state 
of stress, and stress gradient are neglected, 

Mr. Lee’s assumption of a parabolic stress-strain diagram allows predic- 
tion of theoretical values of strength and deformation of reinforced concrete 
members under short duration static loads in a relatively simple manner, 
However, additional experimental evidence is needed to establish the validity 
of the parabolic stress-strain relationship, particularly for predicting defor- 
mations of beams at loads approaching ultimate strength. 


2. “Flexural Behavior of Reinforced Concrete Beams”—by Stuart H. Barthol- 
omew, M.S. Thesis, University of California, Berkeley, California, 
October 1952, 429-4 
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STRESS-STRAIN CURVES FOR CONCRETE 


A. J. ASHDOWN®-— The results of Mr. Lee’s estimation of concrete 
stresses in a beam would certainly be erroneous if obtained on an ordinary 
reinforced beam, since the interaction of the steel with the concrete is not 
as ordinarily assumed. The curve of stress obtained by J. M. Prentis was 
on a non-grouted prestressed beam and was probably nearer to the actual 
conditions. 

In my opinion, the best way of obtaining concrete stress distribution is to 
carry out tests on beams reinforced with bars anchored at the ends and free 
to move in holes, The beams should have stirrups passing around the holes 
to prevent shear failures. Even so, the curve produced would only represent 
the history of the top fibre. 

It is significant that Mr. Lee’s readings stopped at 84 1/2% of the failing 
load; in this interval a considerable amount of strain would have developed 
with a loss of extreme fibre stress, thereby throwing more stress on the 
layers below the top. 

If the fundamentals of Mr. Lee’s paper are invalidated, the complicated 
formulae produced are unwarranted and frightening to the ordinary engineer, 
who wants something simple. 


3. Civ. Eng. Dept., Imperial College, London, England, 
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DISCUSSION OF IMPULSIVE MOTION OF ELASTO-PLASTIC BEAMS 
PROCEEDINGS—SEPARATE NO, 287 


P. S. SYMONDS!—The authors have made the valuable contribution in this 
paper of showing how the familiar normal mode type of analysis of elastic 
structures can be extended to problems where plastic flow must be considered 
in the transient response to a dynamic load. 

There is no doubt that problems of dynamic loading of beams and frames 
into the plastic range are very difficult ones, so that the assumptions made 
by the authors appear necessary to render the problems tractable to analysis 
of moderate length and complexity. Their basic assumptions concern the 
physical properties of the beam, and can be summed up by saying that the 
moment-curvature relation has the same shape as the stress-strain diagram 
of their Fig. 1; thus it is assumed to consist of two straight lines, one at the 
elastic angle with slope EI, and the other parallel to the curvature axis and 
lying a distance Mo above this axis, When this “capacity moment” M, is 
maintained at a particular section curvatures of arbitrarily large magnitudes 
can occur there, and “plastic hinge” action can take place, 

Even with these assumptions the analytical and numerical work is often 
lengthy, and shorter approximate methods become highly desirable. The au- 
thors’ methods for obtaining an upper bound (p. 13) and a lower bound (p. 15) 
are examples of such approximate methods. 

The writer has been concerned with applications of the so-called “rigid- 
plastic” type of analysis of dynamical problems of structures, in which the 
basic simplification is made that deformations of elastic magnitude are ne- 
glected altogether. This type of analysis should give satisfactory results for 
the final deformations when the work done in plastic deformation is very much 
greater than the maximum elastic energy that could be stored in the beam. 
On this basis solutions for a variety of dynamical problems of structures have 
been found, in most cases with relative ease: problems of impact, of impul- 
sive motion (initial velocities) and of various types of specified force-time 
curves have been solved for beams; a problem of a circular ring under con- 
centrated side load has been treated; and certain problems of transverse 
bending of circular plates have been worked out, 

The author’s upper bound solution is obtained on the assumption that defor- 
mation occurs only at a hinge at the mid-point of the beam, This does not 
give the correct solution according to the rigid-plastic hypothesis, Their 
upper bound solution would be the correct rigid-plastic solution only for a 
beam strengthened in the interior of each half so that no deformation could 
occur except at the mid-point where the capacity moment is M,. The correct 
rigid-plastic analysis for the uniform beam takes proper account of the plas- 
tic flow which occurs at sections other than the middle one. This analysis 
has been worked out recently by Cotter(1), The results for the central angle 
for the cases worked out by the authors are shown in Table I. It is seen that 
the final value of 6* from the correct rigid-plastic analysis is considerably 
less than the final value of 6 (from the authors’ elasto-plastic analysis) and 
of 6y (given by the authors’ upper bound solution). 


1. Associate Prof. of Eng., Brown University 
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Table I 


Final Angle et Mid-Point of Free-Free Beam 


NRW 
Vwi 


on 


Correct Plestic-Rigid Solution (Ref.[1] ) 
Upper Bound Solution 
Elastic-Plestic Solution Bleich and Salvadori 


Lower Bound Solvtion 


Similarly, in their elasto-plastic solution for the uniform beam problem 
which the authors have worked out, they have taken proper account of plastic 
flow occurring at the mid-point of the beam, but have not considered possible 
violations of the yield condition at other points, at least in the paper. Again, 
the final deformations they compute would be correct only (in general) for a 
beam strengthened in the interior of each half-beam so that plastic flow could 
occur only at the mid-section, and hence their values of the central angle are 
in general too high, Although the authors state (page 9) that the solution for 
the elasto-plastic motion in terms of normal mode vibrations could be car- 
ried through for cases other than that where a single stationary hinge occurs, 
it would appear that the practical numerical difficulties would be very great. 
This perhaps explains why they have not included in the paper a discussion of 
methods for examining all cross sections of the beam to ensure that the ca- 
pacity moment is not exceeded, 

The authors have compared their single-hinge elasto-plastic solution with 
their single-hinge “rigid-plastic” solution, and find that for the value v,=11.1 
the final angle values differ by about 20 per cent, while for the high value 
Vo = 67.4 the final results are identical, Since both their solutions assume 
plastic flow localized at a single hinge, the conclusions stated in their final 
paragraph seem untenable. In the first place, they have not compared a true 
rigid-plastic analysis with their elasto-plastic solution. More importantly, 
since their single-hinge “rigid-plastic” solution is certainly wrong, the ob- 
served agreement between the final angle computed by their upper bound 
method and that given by their elasto-plastic analysis shows that the elasto- 


429-8 


|_| 
Vo = 11.1 Vo = 67.4 | 
0 
205 0 236 
42 93 
.20 1.32 
on 095 1.46 
1.0% 1.53 
035 1.07 1.56 
+0 1.08 1.58 
1.08 1.58 | 
| 
o* - 
{ 
Sy - 
Op - 
i 


plastic solution likewise is erroneous. The agreement between their two 
solutions for large initial velocities is to be expected, since both apply to fic- 
titious beams strengthened in such a way that plastic flow occurs only at the 
mid-section. The point is that this comparison gives no indication of the cor- 
rectness of either solution, and while their elasto-plastic solution gives cor- 
rect results at low velocities, both this and their upper bound solution are 
erroneous at large velocities. 

The true rigid-plastic solution obtained by Cotter(1) can be expected to be 
correct at sufficiently high initial velocities. According to the criterion pro- 
posed by Lee and the writer(2) results obtained by the rigid-plastic type of 
analysis should be good approximations if the work done in plastic deforma- 
tion greatly exceeds the maximum possible elastic strain energy that could be 
stored in the beam. The ratio of these two energies is about 1 in the authors’ 
case where Vv, = 11.1, and about 36 in the case of v, = 67.4. Since valid 
elasto-plastic solutions for a uniform beam problem have not yet been ob- 
tained one cannot yet say how big this ratio must be in order for the rigid- 
plastic analysis to give satisfactory results. Comparisons have been made 
for simple mass-spring systems, however (3), (4), That given in(4) indicates 
that when the energy ratio is higher than 20 the rigid-plastic type of analysis 
gives results accurate to within 5 per cent in an initial velocity problem, It 
seems fairly likely from this result that the rigid-plastic solution obtained 
by Cotter for the initial velocity 67.4 is reasonably accurate, 

Reference(4) deals with perhaps the simplest impulsive motion problem of 
a uniform beam, namely that of a simply supported span with an initial veloc- 


ity proportional to sin | x, { being the span length, It is easy to verify di- 


rectly in this problem that a “single-hinge” type of elasto-plastic analysis is 
erroneous even for fairly small initial velocities. Even in this problem, how- 
ever, it seems to be very difficult to obtain the correct elasto-plastic solution 
for arbitrary initial velocities, It would be most valuable to have the benefit 
of the authors’ experience in dealing with such problems, 
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DISCUSSION OF DESIGNING ALUMINUM ALLOY MEMBERS 
FOR COMBINED END LOAD AND BENDING 
PROC!!EDINGS—SEPARATE NO. 200 


R. R. McCALLEY, Jr.(1) J, M, ASCE.—Unsymmetrical cases occur suffi- 
ciently often in practice that a generalized form of Eq. (6) may be of interest. 
Eq. (6) applies to the case in which the constant bending moment lies in one 
of the principal planes of the cross section. Consider the more general case 
in which the x and y axes do not coincide with the principal axes, For a 
beam subject to a constant bending moment about the x axis, the critical 
moment, Mery, can be expressed as: (2 


1, (KL)? GJ 


where I, is the major moment of inertia, Ig is the minor moment of inertia, 
and Ix is the moment of inertia about the x axis. 
However, the generalized definition of e is now given by: 


Fi * Fu (b) 


where y; is the y coordinate of a geometric point designated the “twist center,” 
and y, is the y coordinate of the shear center. The “twist center” is the 

point about which a beam in pure torsion rotates. In symmetrical cases both 
the “twist center” and the shear center coincide with the centroid of the cross 
section, The “twist center” coordinates (x;,x;) may be computed from the 
formulas; (2 


1 2 2 


2 2 


where P, is the radius vector from some convenient reference point (xg,ya) 
to a point (x, y) on the cross section and the integrations are taken over the 
whole area, A, of the cross section, 

The critical stress f, may then be determined from the formula: 


1, Engineer, Applied Physics Laboratory, Johns Hopkins University, Silver 
Spring, Maryland. 


2. “Contributions to the Theory of Combined Flexure and Torsion” by 
R. B. McCalley, Jr., Thesis, Cornell University, June 1952, 
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The x and y coordinates of the extreme fiber must be chosen by inspection 
for the worst case, In many instances, it may be convenient to eliminate 
the product term I; Ig in the foregoing equations by means of the identity: 


2 
1,1, * 1, if) 
In the event that the x axis coincides with the major axis of the cross sec- 
tion, L, = Ij, ly = Iz, Ixy = O, and Eq. (a) reduces to Eq. (6). Goodier (3) has 
previously called attention to the fact that the = sign before the radical term 
indicates that the beam has one critical moment when bent downward but an- 
other when bent upward. The quantity e is equal to one-half of Timoshenko’s 
B,.(4) Notice that the expression given by Timoshenko is referred to the 
principal axes of the section. 


3. Reference 9, page 10. 


4. “Theory of Bending, Torsion, and Buckling of Thin-Walled Members of 
Open Cross Section” by S. P. Timoshenko, Journal of the Franklin Institute, 
May 1945, page 356. 
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